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We find that several thresholds can contribute to the enhancements of the newly observed heavy 
pentaquark candidates P(l*'(4380) and P^(4450) via the anomalous triangle singularity (ATS) tran¬ 
sitions in the specific kinematics of Aj, —>■ J/^|lK~p. Apart from the observed two peaks we find 
that another peaks around 4.5 GeV can also be produced by the ATS. We also show that the 
can be produced at leading order in A^ decay. This process is different from the triangle diagram 
and its threshold enhancement only appears as CUSP effects if there is no pole structure or the 
ATS involved. The threshold interaction associated with the presence of the ATS turns out to be a 
general phenomenon and plays a crucial role in the understanding of candidates for exotic states. 

PACS numbers: 12.39.Mk, 14.20.Pt, 13.30.Eg 


I. INTRODUCTION 

The two states P+^380) and P+(4450) observed in the invariant mass spectrum of J/ifp in Af, —>• J/'ipK~p by the 
LHCb Collaboration [l| has immediately attracted a lot of attention from the whole community since they could be the 
long-searching-for pentaquark states in the heavy flavor sector. Recent theoretical studies can be found in Refs. [H- 
Their masses are 4380 ± 8 ± 29 MeV and 4449.8 ± 1.7± 2.5 MeV, respectively, and their widths are 205 ± 18 ± 86 MeV 
and 39 ± 5 ± 19 MeV, respectively. Their preferred are either 3/2“ and 5/2+, or 3/2+ and hj^T, respectively, 
based on the detailed experimental analysis. The possible existence of multiquark states has always been regarded 
as a natural consequence of QCD. Although the conventional quark model has made tremendous successes in the 
description of the hadron spectroscopy, it also raised questions on why and how those multiquark states kept out of 
our sight for such a long time. The LHCb results and the recent results from Belle Q and BESIII @ certainly make 
a big step forward to our understanding the multiquark system. But at the same time, they also give us chances to 
ask more questions. 

In this Letter we will examine the role of the non-perturbative anomalous triangle singularity (ATS) in the decay 
of Ab —>■ Jl'f}K~p. As recently pointed out in Ref. the pronounced narrow threshold states in the elastic channels 
should indicate non-perturbative rescatterings which will eventually generate pole structures after sum over all the 
loops to infinity. It was also stressed that if the ATS is present the threshold peak will be enhanced and mix with 
the dynamic pole structure in the inelastic channel. Therefore, in order to understand the nature of the threshold 
enhancements a careful analysis of the triangle process should be necessary. 

There are several interesting features arising from the decay channel as the data have shown. For instance, there 
are clear structures for A* resonances in the lower end of the K~p invariant mass spectrum while the higher mass 
region appear to smooth out. In contrast, the observed P+(4380) and P+(4450), in particular, the P+(4450), are 
clear structures above the phase space. Note that some interferences occur over a rather broad mass range above 
the P+(4450) in the J/ifp invariant mass spectrum (see e.g. Fig. 2 of Ref. i). This feature suggests that there 
are intermediate processes in A;, —>■ J/'tpK~p which will give rise to final state interactions. Our motivation is to 
investigate whether the ATS will accumulate the events at the two peak positions. 

As follows, in the next Section we will first analyze the transition mechanisms for Af, —>■ J/'!/Ar“p and then investigate 
the loop diagrams where the ATS can be present. A summary will be given in the last Section. 
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FIG. 1: Relevant processes for the production of P^(4380) and P^(4450) in At, —>■ J/'ijjK p. 


II. PRODUCTION MECHANISMS FOR THRESHOLD STATES 


The transition of A{, —>■ J/'ij}K~p is dominated by the flavor changing weak decay of 6 —?► c + cs via the V — A 
current. This leads to an intuitive expectation that the u and d quark in should be a spectator, hence their isospin 
quantum number / = 0 should be conserved. As a consequence, if the ud pair is to be combined with a c quark to form 
a charmed baryon, it will favor A^ baryons instead of Ec ones. In Refs. the authors propose that P+ (4380) 

and P+(4450) are molecular states of Ec(2455)ij*(2007) and Sc(2520)Z)*(2007). But it was not discussed why the 
production of Ec*^ in the A^ decay should be significant. Also, in Ref. [I^ the authors refer one of the pentaquark 
candidates to a molecular state of Ec(2455)i5*(2007) with = 3/2“. In Ref. [l^ the authors proposed to produce 
the Pi*" (4450) via the J/ipp rescattering to D^*'^ which, however, should be higher order terms. In a followed-up 
paper [ij the authors investigate other isospin 0 channel recoiled by J/'0. 

Actually, the production of the E), is not necessarily to be suppressed at leading order. Since the b quark decay 
is a short-distance process, the final c, c and s quark will be created locally. They can easily split the ud quark and 
rearrange them with the created uu to form, e.g. K~. Then the rescattering of D^*'> into J/ipP can occur. 

This process is illustrated in Fig. [T](c) with (a) and (b) the same mechanism analyzed in Ref. [I| as a contrast. 

For Fig. [11(b), a decomposition of the isospin gives the relative decay amplitudes 

{Y,KD\H^\Kb)ib) 

= ^(A+A0P--A+A-P0) , (1) 

where Yc denotes charmed baryon Ac*^ or Ec*^ We denote the possible Ac*^ and by the ground state symbols 
and simply write the amplitudes by the final state particles. The flavor wavefunction {ud — du)b/\/2 for Af, is implied. 

Similarly, the decomposition for Fig. [T](c) can be written as 


(Y;AP|i7^|A6)(,) 

1 


2y/2 


-Yt^K-D' 


\YtK^D- - \YtK- 


+ Y°K°D° 


-AtK-D° - -A+K°D~ 
2 2 


( 2 ) 


where Y'^^K~D~ and Y'^K~D^ would allow the formation of pentaquark states via the Y^'^D^*'^ rescatterings as 
studied in Refs. (l(Tl - [r^ . From Eqs. (HJ and (I2|) we can see that Fig. |T](b) dominantly contribute to the production of 
7=1/2 pentaquark system recoiled by K~. 

If the production of the pentaquark candidates P+(4380) and P+(4450) is indeed via the Ec*^75^*^ fTable HV)) or 
A(*)/)(*) (Table [n]) interactions, then the relevant thresholds will play non-trivial roles in the coupled-channel decays 
into e.g. J/'4’P- Considering the major rescattering processes arising from Fig. [1] we can reexpress the transitions 
as Fig. [2] where three type of rescatterings will contribute. Figure [5] (a) and (b) are triangle diagrams while (c) 
is two-point loop interactions. We note that the creation of E^*^ via Fig. [T] (c) should be driven by short-distance 
interactions. This process will favor the formation of (4380) and (4450) if they are states generated by Ei*^7)C) 
interactions. 












3 



FIG. 2: The loop diagrams as a consequence of Fig. [T] where the ATS and kinematic CUSP can be recognized. 


The interesting property of Fig. [5] (a) and (b) is that given the masses of the involved states located within certain 
ranges it will allow the internal states to be on-shell simultaneously. This is different from the kinematic CUSP 
effects which only recognizes the on-shell condition for two internal particles and contributions from such a branch 
point is subleading compared to the ATS [1^ [l^. Therefore, we do not expect that such CUSP effects produce 
narrow and strongly enhanced structures in the invariant mass spectrum. In contrast, when the ATS coirdition is 
satisfied, the siirgularity behavior of the integral will produce strong enhancements at the singular points of which the 
effects can be measured in the experiment. In particular, the singular points will mostly locate in the vicinity of the 
two-body thresholds but not necessarily to be exactly at the thresholds. It should be realized that the positions of 
the singularity will not change even when higher partial waves contribute at the interaction vertices.^ The reason is 
because the singular term will always be kept in the decomposition of the integrand in the Feynman parametrization. 
In other words, even though the contribution from the singular term relative to other contributions might be small, 
its enhancement at the singular point may not be negligible^. Nevertheless, in the case of A;, —>■ there are 

several thresholds close to each other. Even a small singularity enhancement can build up and produce measurable 
effects. 

Since quite a lot of thresholds can appear in the decays of Fig. [2] and we are still lack of information about the 
vertex couplings, we only consider low partial waves and thresholds which are close to the masses of interest and we 
discuss separately the properties of those three types of loops in Fig. [21 

Figure 121 (a) is a consequence of Fig. [1] (a) where the rescattering between Xcj and an exchanged proton from the 
decay A* —>• pK~ is considered. Note that the mass thresholds for p + Xcj {J = 0, I, 2) are close to the peak 
masses for P+(4380) and P+(4450) as listed in Table [H Also, the S-wave scatterings of pXc 2 J/ipP ^ can access 
the quantum numbers of 3/2“'' and 5/2’'' for the threshold enhancement. The Xci and p scattering can access the 
quantum numbers of 1/2''" and 3/2“*'. The XcoP can reach 1/2“ and 3/2“ via a P wave interaction. It is interesting to 
notice that the significant enhancement to the XcJP (with J = 0, 1, 2) via the ATS would prefer that the mass of A* 
within the mass regions 1.92 ~ 2.20 GeV, 1.89 ~ 2.11 GeV, 1.83 ~ 2.06 GeV, respectively. From Fig. 2 (a) of Ref. [l|, 
it shows that the cross section for K~p is smooth but non-zero. Note that as long as the kinematics approaching the 
ATS condition, all the cross sections will contribute to the threshold singularity. 

In Fig. [3] we show the structures in the invariant mass of J/4’P via the triangle diagram of Fig. |2] (a). As discussed 
before, since XciP and Xc 2 P can access the possible quantum numbers via the S wave, we only consider loops of Xci 
and Xc 2 at this moment. On the other hand, since the branching ratio of —>■ K'^Xc 2 is one order of magnitude 


^ The only difference is that the higher partial waves will make the signal from the ATS less pronounced. 

^ The detailed discussion about the ATS and their manifestations in physical processes can be found in R ef. |15|| and there are cases that 
the ATS involving higher partial wave interactions can still produce significant threshold enhancements [ld42lll . 

® Although the production of XcO 3-ud Xc 2 is suppressed in B decay or any hadron with only one heavy b quark without charm or anti¬ 
charm quark as discussed in Ref. 1221 . we can still estimate how large the next leading order contributions are or what is the behavior 
after considering it. 
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TABLE I: The x^jP thresholds which can be enhanced by the ATS via Fig. [3 (a). 


Threshold masses [GeV] 

Xco(lP) 0+ 

Xci(lP) 1+ 

Xc2(lP) 2+ 

pl/2+ 

4.353 

4.449 

4.494 



FIG. 3: The invariant mass distribution of J/'i/’P given by the triangle diagram of Fig. [5] (a). The vertical dashed lines indicate 
the masses of PJ*’(4380) and P^(4450), respectively. Here the contribution of Xc 2 is estimated as the order 10~^ of that from 
Xcl- 


smaller than that of 5+ —>■ K'^Xci, we expect that it is also the case in the A;, decay. Then, the upper limit of the 
contributions from Xc 2 can be estimated by requiring BR{Ki, —>■ K~Xc 2 p)/BR{Ai, —>■ K~Xcip) ~ 10“^. As shown in 
Fig. [21 the ATS can produce significant threshold enhancement of XciP while the effect from the Xc 2 P loop is strongly 
suppressed. Note that the singularity from the XciP loop exactly locates at the mass position of the observed Pc (4450). 
The pole trajectories of the XciP loop and the properties of the poles on different Riemann Sheet were first considered 

in Ref. [13]. 

It is possible that the intermediate cs in Fig. [T] (b) can form intermediate d[*j states which can decay into 1)1*1 
and K~. The intermediate Z)l*l meson will then scatter the Ac into J/ipp. This process is illustrated in Fig. [2] 
(b). The accessible thresholds are listed in Table [TTl Although one can see from Table [IT] that in an S wave none of 
the thresholds matches the experimental measured masses and favored quantum numbers simultaneously. However, 
among these ‘Ac’s, the Ac(2595)+i)sj(2860) loop with i)sj(2860) decay to DK can produce the singularity at the 
P+(4450) mass position as shown by the black point in Fig.|2|(b). As a comparison we also show the pole trajectory 
with i)sj(2860) decay to D*K in Fig. Uj (a). Note that the ATS only works in a very limited kinematic region. It 
implies that if the kinematics deviate from the ATS condition one should not expect any significant enhancement at 
the corresponding threshold mass region. This will provide a possibility for experimentalists to pin down the nature 
of some threshold states. Namely, if they are not caused by kinematic effects the enhancements will still appear in 
other processes where the ATS condition does not hold. 

Although Ac(2595)i)sj(2860) can access the expected quantum number in a P wave which will be suppressed by 
the centrifugal barrier, we will show later that the singular points near threshold can still match the observed peak 
positions. Also, as mentioned earlier that the ATS can still possibly produce observable effects when higher partial 
waves are present at the interaction vertices, we then investigate the possible partial waves for Fig. [5] (b) and see how 
the ATS would manifest in the invariant mass spectrum. 

In Table Hill the thresholds for the Ai*^D^*j around the mass of A^ are listed. To fully satisfy the ATS condition, 
the threshold of Ai*'^ D^*J should be close to the mass of A;,. We note that for most of the Ac*^ thresholds, they 


have already allow a sizeable enhancement near the Ai*'^D^*'^ thresholds. 

In Fig. 121(b), the charmed baryon and anticharmed meson are chosen as Ac(2595)/Ac(2650) and D/D* , respectively, 
which means that we take into account the contributions from four diagrams. Furthermore, to match the qualitative 
feature of the experimental observations, we estimate the rescattering amplitudes by assuming relatively smaller 
branching fractions of A^ Ac(2625) and D^*] — D*K, compared with those of Af, —>• P^jAc(2595) and 

D^*j — >■ DK^ respectively. The numerical results of J/'ipP invariant mass distributions given by the diagrams of Fig. | 2 ] 
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FIG. 4: The trajectories of the two singularities of Fig[2] (b) in terms of the mass of dI*J with tuac = 2.595 GeV. The 
emitted charmed mesons from D^*J are D* (a) and D (b), respectively. The masses of d[*J decrease along the arrows. In 
order to distinguish the pinch singularities, a small width, 10 MeV, is assigned to d[*J . The black point in (b) corresponds to 
the Dsj(2860). The red thick solid lines are the possible mass regions which can produce the ATS enhancement. The vertical 
dot-dashed lines in both plots indicate the mass position of the (4450). The corresponding trajectory of Fig. (a) has been 
investigated in Ref. [^ . 


TABLE II: Thresholds accessible in the invariant mass spectrum of J/ipp. The two numbers in the square bracket have beyond 
the allowed phase space for J/ipp. 


Threshold masses [GeV] 

Ac(2286) 1/2+ 

Ac(2595) 1/2“ 

Ac(2625) 3/2“ 

Ac (2880) 5/2+ 

U(1865) 0" 

4.151 

4.457 

4.493 

4.746 

U*(2007) 1' 

4.293 

4.599 

4.635 

4.888 

L>i(2420) 1+ 

4.706 

5.015 

5.045 

[5.300] 

792(2460) 2+ 

4.746 

5.055 

5.085 

[5.340] 


(b) are displayed in Fig. There will be many options for D''*j, we therefore display the results for several masses 


and widths of D^*j. Since the intermediate state D^*j in Fig. [2Ib) could be broad, to account for the width effect 
we adopt a Breit-Wigner-type propagator [q^ — -|- im^(.)F = (.)]“^ in the calculation of the loop integrals. The 

^sj ^sj ^aj 

complex mass of an intermediate state will move the ATS from the physical boundary by a small distance. If the 
width is not extremely large, the ATS will lie close to the physical boundary, and the scattering amplitude can still 
feel the influence of the singularity. In Fig. [51 when the mass of D^*j is taken around 2.86 GeV, one notices that 
four peaks arise in the J/ipp distributions, which stay in the vicinities of Ac(2595)Z), Ac(2625)Zi, Ac(2595)l)* and 
Ac(2625)I)* threshold, respectively. The narrow peak around the Ac(2595)Z) threshold can match the structure of 
Pc(4450) as observed in experiment. When the widths of vary from 0 up to 150 MeV, the main feature of the 
threshold enhancements still hold though the rates will decrease. 

For Fig. [2] (c), the kinematic effects will be just CUSP structures in the J/ipp invariant mass spectrum. We do 


not discuss the dynamic consequences if the intermediate D^*'^ (Table HVl) may have strong couplings, then they 
may generate dynamic poles near threshold after proper summation over the bubble loops. Instead, we only show 
the kinematic CUSP for which has turned out not to lead to pronounced structures as recently studied in Ref. 

In another word, the observed pronounced peaks may either be produced by possible pole structures or the ATS 
mechanisms. 

In order to try to distinguish the behavior from a pole structure and the ATS, we generate the invariant mass 
spectra of J/ipP in different K~p invariant mass regions the same as Fig. 8 of Ref. [I[ (see Fig. [6]). Since the ATS 
contributions will vary in terms of different kinematics, such a quantity will be able to distinguish the behavior of the 
pole and the ATS mechanism. The results for those three loop processes of Fig. |2] (a), (b) and (c) are illustrated by 
the solid, dashed and dotted lines which are different from the symmetric Breit-Wigner lineshape. In particular, the 
structures created by the CUSP effects appear to be negligible. Note that the kinematics of higher invariant mass of 
K~p is favored by Fig. [5] (a). It means even small couplings for XcJP scattering are introduced, the ATS enhancement 
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TABLE III: Thresholds for the d[*J Ai*'>. 


Threshold values in the square bracket are above the mass of A6(5619). 


Threshold masses [GeV] 

Ac(2286) 1/2+ 

A.(2595) 1/2- 

A.(2625) 3/2" 

A42880) 5/2+ 

L>41968) 0" 

4.254 

4.563 

4.593 

4.848 

U:(2112) 1- 

4.398 

4.707 

4.737 

4.994 

5.0(2317) 0+ 

4.585 

4.912 

4.942 

5.197 

5.1 (2460) 1+ 

4.728 

5.055 

5.085 

5.340 

5.1 (2536) 1+ 

4.822 

5.131 

5.161 

5.416 

5.2(2573) ?■ 

4.859 

5.168 

5.198 

5.453 

5.1 (2700) 1" 

4.986 

5.295 

5.325 

5.580 

5.j(2860) ?• 

5.146 

5.455 

5.485 

[5.740] 

5.j(3040) ?• 

5.331 

[5.636] 

[5.672] 

[5.926] 


TABLE IV: The thresholds accessible in the invariant mass spectrum of J/ipP- The two numbers in the square bracket 

have beyond the allowed phase space for J/'4’P- 


Threshold masses [GeV] 

Eo(2455) 1/2+ 

Ec(2520) 3/2+ 

E<,(2800) ?■ 

5(1865) 0" 

4.321 

4.385 

4.668 

5* (2007) 1" 

4.463 

4.527 

4.810 

5i(2420) 1+ 

4.875 

4.939 

[5.222] 

52(2460) 2+ 

4.917 

4.981 

[5.264] 


can still be observable. 

Since we introduce several mechanisms to generate the kinematic singularities, it is also necessary to discuss their 
similar and different characteristics. Both XciP and Ac(2595)l) thresholds are very close to the mass of PT(4450). 
Note that XciP can scatter into J/'i/'P via multi-gluon exchange process, while A*Z) can scatter into J/'f/’p via quark 
interchange process which is a rearrangement of the quark flavors. The P-wave scattering between Xci and proton will 
imply that the quantum numbers of the J/tpp system can be = (l/2“, 3/2“, 5/2“), but the P-wave scattering 
between Ac(2595) and D will imply that the quantum numbers can be = (l/2“,3/2“). Some of these quantum 
numbers are compatible with the experimental fitting results. If one hopes the rescattering mechanism can favor the 
= 5/2+ assignment of P+(4450), the P-wave rescattering of XciP or Ac(2595)P will be required. Usually the 
higher partial wave rescattering amplitudes will be suppressed to some extent. However, the analytic properties of the 
kinematic singularities will mainly depend on the kinematics of the loop integrals, which will not be affected too much 
by the coupling forms of the vertices. To clarify whether the XcJP rescatterings or the rescatterings would 

be dominant in producing the resonance-like peaks, we suggest that A° — K~hcP would be a promising channel. 
Namely, the transition XcJP ^cP would require the heavy quark spin flip, thus, will be suppressed. In contrast, in 
A*D^*)^hcP such a spin flip does not necessarily occur. Furthermore, since the hcP threshold is much larger than the 
lower thresholds, some higher resonance-like peaks induced by the rescatterings, similar to P+(4450), can be expected 
in the hcP invariant mass distributions. 


III. SUMMARY 

In this work we analyze the role played by anomalous triangle singularity in A^ —>■ J/il)K~p as a possible contribution 
to the observed hidden-charm pentaquark candidates P+(4380) and P+(4450). We first show that the D^*'> can 
be produced at leading order which seems to be overlooked in the literature. We then demonstrate that the ATS 
can generate threshold enhancements which can mimic the experimental observations. We have to admit that more 
detailed model construction is needed in order to determine better the ATS behavior and strengths. But we emphasize 
that the kinematics of Af, —^ J/'0A“p are in favor of the ATS mechanism when intermediate states are involved. Since 
all the triangle amplitudes can build up the contributions from the ATS near the mass region that we are interested 
in, they seem to be able to produce structures similar to what observed in experiment. Meanwhile, we hnd that the 
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^7/iAp [GeV] 

FIG. 5: Invariant mass distributions of J/ipp given by the triangle diagram of Fig. [3 (b). The width of the intermediate d''*j 
are taken as (a) 150 MeV, (b) 50 MeV and (c) 0 MeV, respectively, as a demonstration of the width effects. The masses of the 
D^*j are 2860 MeV, 2700 MeV, 2573 MeV and 3040 MeV corresponding to D^*j states listed in Table [ITT] The vertical dashed 
lines indicate the masses of P,^(4380) and ^,^(4450), respectively. 


Ac(2625)£>(1865) threshold can give rise to a third peak around 4.5 GeV. 

We also discuss the kinematic feature of which can produce kinematic CUSP effects. In this sense, the 

observed pronounced structures can be either produced by pole structures or the ATS mechanism. The kinematic 
dependence of the ATS is investigated by looking at the ATS cross sections in different K~p invariant mass regions. 
We also show that the ATS mechanism should have rather strong dependence of the kinematics by looking at the 
J/ipP invariant mass distributions at different energy cuts for K~p. 

In brief, since the ATS could play a role in the production of threshold states, it may mix with the pole if the 
singular threshold is close to the pole position. Such a mixing and interference might generate complicated structures 
such as to distort the lineshape or produce narrower enhancements that behaves very differently from expected pole 
contributions. In order to better understand the nature of these newly observed pentaquark candidates, a combined 
study of the ATS and dynamically generated pole structures should be necessary. 
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FIG. 6: Invariant mass distribution of J/ijjp with different K~p momentum cuts which are same as Fig. 8 of Ref. [I|], i.e. (a) 
rriKp < 1-55 GeV,(b) 1.55 GeV < niKp < 1-07 GeV, (c) 1.07 GeV < niKp < 12.0 GeV, (d) mxp > 2.0 GeV. The vertical 
dashed lines indicate the masses of Pj''(4380) and Pt,’*"(4450), respectively. 
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